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Gravitational Depolarization of Ultracold Neutrons: Comparison with Data 
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We compare the expected effects of so-called gravitationally enhanced depolarization of ultracold 
neutrons to measurements carried out in a spin-precession chamber exposed to a variety of verti¬ 
cal magnetic-field gradients. In particular, we have investigated the dependence upon these field 
gradients of spin depolarization rates and also of shifts in the measured neutron Larmor precession 
frequency. We find excellent qualitative agreement, with gravitationally enhanced depolarization 
accounting for several previously unexplained features in the data. 

PACS numbers: 13.40.Em, 07.55.Ge, ll.30.Er, 14.20.Dh 


I. INTRODUCTION 

Ultracold neutrons (UGN) are neutrons of extremely 
low energy, typically 200 neV or less, which can be stored 
in material bottles and which are routinely used in exper¬ 
iments such as the ongoing search for the neutron electric 
dipole moment (nEDM). Gollisions with the containing 
walls are elastic, so the UGN never thermalize. Being 
of such low energy, they “sag” under gravity, and rather 
than being distributed uniformly throughout their stor- 
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age vessel their density decreases with increasing height, 
with each specific energy group having its own center of 
mass. In the presence of a vertical magnetic-field gra¬ 
dient, the average magnetic field sampled by the neu¬ 
trons will therefore depend upon the neutron energy. The 
implications of this stratification have been discussed in 
earlier work BH], but, in summary, it results in a rel¬ 
ative dephasing of the neutrons in different energy bins, 
which then alters the measured Larmor spin-precession 
frequency. This phenomenon is referred to as gravitation¬ 
ally enhanced depolarization, in contrast to the intrinsic 
depolarization that takes place within each energy bin as 
a result of the neutrons sampling different fields as they 
move around the storage volume. A key distinction is 
the asymmetric nature of the gravitationally induced de¬ 
phasing, as shown in Fig. 3 of [T]) with the lowest-energy 
neutrons playing a particularly crucial role. The resulting 
nonlinearities in frequency response as a function of ap¬ 
plied magnetic field gradients represent potential sources 
of systematic uncertainty in precision experiments such 
as nEDM searches [Mg. Since such experiments pro¬ 
vide tight constraints on physics beyond the Standard 
Model, with consequent implications for particle theory 
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and cosmology, a full understanding of the phenomenon 
is essential. 

In this article, we compare our experimentally mea¬ 
sured results, both in terms of frequency shifts and of 
depolarization rates, with those anticipated from theo¬ 
retical calculations. We begin in Section [iT] with a dis¬ 
cussion of the spectrum of UCN within our storage cell; 
this underlies the subsequent calculations of the gravita¬ 
tionally enhanced depolarization. We give an overview of 
the calculations themselves in Section Hill In Section ITVl 
we discuss the basic intrinsic-depolarization mechanisms, 
which are revealed to make only a minor contribution to 
the frequency shifts. We then present, in Section a 
direct comparison of the anticipated and measured po¬ 
larization a remaining after 180 s of storage in a range 
of applied S-field gradients. In Section |VI| and |VII| we 
consider the frequency shifts that arise from this phe¬ 
nomenon, before finally discussing in Section |VIII| the 
possible implications for nEDM experiments, including 
the current world limit in particular. 

The measurements described in this article form part 
of a program of work aimed at an accurate determina¬ 
tion of the nEDM, currently being carried out at the new 
high-intensity UCN source [7] based at the Paul Scherrer 
Institute (PSI). The experimental apparatus and proce¬ 
dures are described in substantial detail in |S]. The ap¬ 
paratus is based upon that used [H] in an earlier nEDM 
measurement at the Institut Laue-Langevin (ILL) |3] , but 
substantially upgraded with the incorporation, in partic¬ 
ular, of an array of Cs magnetometers m, a system for 
the simultaneous detection of both neutron spin states 
m, and a set of active compensation coils that provide 
dynamic shielding of external magnetic fields |12| . 

The 1 /rT magnetic holding field Bq within the EDM 
spectrometer is primarily vertical, so Bq B^, although 
there are small transverse components B^, By present at 
the ~ few nT level. We define 

+ ( 1 ) 


where B^ = B^ + By. 

In order to compensate for changes in Bq, a cohabiting 
atomic mercury magnetometer m is used to make pre¬ 
cise real-time measurements of the volume-averaged field 
within the UCN storage cell. Under an applied vertical 
magnetic-field gradient dBz/dz, the measured ratio R of 
neutron to mercury precession frequencies undergoes a 
relative change of, to first order. 
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where Ah is the difference between the centers of mass of 
the populations of (thermal) mercury atoms and (ultra¬ 
cold) neutrons. Precise measurements of this frequency- 
ratio dependence are the subject of [S]. As we shall 
see, gravitationally enhanced depolarization can impose 
a substantial nonlinearity in this relationship: indeed, we 
are unaware of any other mechanism that can do so to 
the extent required to match our observations. 


II. INPUT SPECTRA 

The extent of gravitational depolarization clearly de¬ 
pends heavily upon the spectrum of stored UCN. We have 
recently carried out a series of measurements using a spin- 
echo technique m , from which we were able to derive the 
distribution of energies of UCN remaining after 220 s of 
storage in our apparatus. The resulting fitted spectrum 
is parameterized by 

p{E) = A- -^-^, (3) 

1 -l- g aeq 1 -|- e 

where A is an arbitrary normalization, Eq = 7.7 neV, 
AEq = IneV, El = 28.7neV, and AEi = 6.25 neV. The 
form of this parameterization is based on a very general 
distribution n{E)dE oc E^/'^dE from the low-energy tail 
of a Maxwell-Boltzmann distribution, allowing for low- 
and high-energy cut-offs. 

The spin-echo technique is particularly sensitive to the 
presence of low-energy UCN, but once the neutrons start 
to populate the bottle more or less uniformly it becomes 
increasingly difficult to distinguish between different en¬ 
ergies. This is clear from Figs. 2(a) and 2(b) in |14| . 
where the low-energy tails are fitted well but the high- 
energy region produces less reliable results. Furthermore, 
the spin-echo measurements were carried out at a stor¬ 
age time of 220 s, whereas the polarization and frequency- 
ratio measurements used in the current analysis were car¬ 
ried out at a storage time of 180 s. On both counts, there¬ 
fore, we should not be surprised if the actual spectrum 
were to be somewhat firmer than that arising from the 
spin-echo measurement. 

We have also used the package MCUCN m to carry 
out a detailed simulation of the UCN within our appara¬ 
tus, which yields an alternative estimate of the spectrum 
after 180 s of storage. The simulation is based upon very 
detailed modeling of the PSI UCN source, beamline, and 
guides, as well as of the nEDM storage vessel. The latter 
consists of aluminum electrodes coated with diamond¬ 
like carbon, which form the floor and roof, and between 
them an insulating cylindrical polystyrene ring coated 
with deuterated polystyrene to provide radial contain¬ 
ment. The simulation accounts for losses during stor¬ 
age both from /3 decay and as a result of wall collisions 
m, with the “loss-per-bounce” factor f = WjV (where 
W, V are the imaginary and real parts, respectively, of 
the Fermi potential) set to a common value of 3 x 10“^ 
for the electrodes and for the insulator walls. In fact, 
although V is well known, W is difficult to determine. 
Losses are likely to be dominated by hydrogen that has 
diffused into the containing surfaces, and can - because 
it has an extremely high incoherent-scattering cross sec¬ 
tion - substantially influence loss rates without signifi¬ 
cantly altering the surface potential V. Using an average 
value of / for all of the containing walls appears a reason¬ 
able approach, and the number here arises from earlier 
simulation-based studies HZ] that were tuned to match 
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the observed numbers of neutrons stored as a function of 
time. 

Any damage to the coatings on either the electrodes 
or the insulating walls would result in an area of re¬ 
duced Fermi potential that would preferentially deplete 
the higher-energy neutrons. The same is true of small 
gaps, which may not be completely accurately modeled in 
the simulation. The actual spectrum, therefore, is likely 
to be somewhat softer than the simulation would suggest. 

The spectra resulting from the Monte Carlo (MC) sim¬ 
ulation and spin-echo (SE) studies are shown in Fig. ^ 
It is useful to refer to UCN energies E in terms of the 
maximum height e = E/{mg) attainable under gravity 
in a trap with no vertical confinement, and to this end 
the abscissa is in units of cm. 



FIG. 1: (Color online) Estimated energy spectrum of stored 
UCN from spin-echo measurement (blue solid line) and from 
simulation (green dashed line). Each is normalized to give 
unit total area under the curve. 

We note here for clarity that the energies are defined to 
be the kinetic energies that the UCN would have at the 
floor of the storage vessel, i.e. at the bottom electrode. 


III. CALCULATIONS 

Calculations of the gravitationally enhanced depolar¬ 
ization effect are relatively straightforward to carry out. 
Phase-space arguments m show that the variation of 
density p with height z of UCN of height-equivalent en¬ 
ergy e is (for z < e) given by 

/ Z\l/2 

= p(0u) (l - -) , (4) 

assuming sufficiently diffuse reflections for the phase 
space to be approximately uniformly populated on a 
timescale short compared to the storage time. (Our 


Monte Carlo simulations confirm that this typically takes 
place within 5 s of closing the UCN shutter if more than 
10% of reflections are diffuse, and more quickly still with 
higher diffusivity and also for the more energetic of the 
neutrons in our spectrum.) The height distribution of 
Eq. I^will then be reflected in the distribution of average 
magnetic fields to which UCN of any particular energy 
will be exposed, from which the distribution (appropriate 
to that UCN energy) of integrated phases acquired after 
180 s of free Larmor precession can be calculated. This 
procedure is carried out for all UCN energies across the 
spectrum, accounting for the relative populations of each 
energy bin. The resulting total array of integrated phases 
(fii is then subject to a Ramsey-type analysis, where, as 
discussed in |T], the net frequency is determined by the 
reference phase 

<^ = 2n7r + tan-i (5) 

V (cos (jii) j 

divided by the Ramsey coherence time (180 s in this case). 
Note that the 2n7r term, which arises because the Ram¬ 
sey technique measures phases modulo 27r, is relatively 
easily accounted for by, for example, monitoring the dis¬ 
crepancy between the reference phase and the mean of 
the array of time-integrated phases, and adding (or sub¬ 
tracting) factors of 27r as appropriate to compensate. The 
nonlinearities in response referred to earlier primarily 
arise when the lowest-energy UCN, which do not reach 
the roof of the storage trap, have an integrated phase 
that differs by more than tt radians from the reference 
phase: they then “wrap around” and appear to enhance 
the high-energy tail of the distribution. We will refer to 
this phenomenon as “Ramsey wrapping”. 

Effects due to intrinsic depolarization, arising from 
both vertical and horizontal field gradients, can also be 
included by appropriate weighting of the distribution of 
phases. We discuss this in some detail in the following 
Section . 


IV. INTRINSIC DEPOLARIZATION 
MECHANISMS 

Detailed calculations of intrinsic depolarization within 
magnetic-field gradients have been carried out elsewhere 
|19H26j . There are four relevant scenarios to consider: 
vertical gradients dB^jdz', horizontal gradients of the 
form dB^/dx] transverse fields Bt and their gradients; 
and wall collisions involving small magnetic impurities. 
We present here some simple and rather intuitive models 
of the depolarization mechanisms, and we discuss possi¬ 
ble implications for the polarization and frequency-ratio 
measurements. Throughout this Section, where calcula¬ 
tions are dependent upon an input spectrum we use that 
derived from the SE measurement. 
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A. Vertical gradients: dBzjdz 

Here we consider UCN confined within a vertical 
magnetic-field gradient dB^/dz. Let the confining trap 
be a cylinder of height H and radius r. Following |TH] , we 
replace H with an “effective height” 'H(e) which is sim¬ 
ply defined as the lesser of e; this accommodates UCN 
with energies too low to reach the roof of the trap. 

The following method is based upon that outlined in 
the derivation of Eq. 68 in m- We shall consider our 
trap to be divided by a horizontal plane into two halves, 
with average field strengths that each differ from the field 
at the center plane by Ai?z = {dB^/dz)!^,/^. Let the av¬ 
erage dwell time for UCN in each half of the trap be tu,. 
Using the standard kinetic-theory result (due to Clau¬ 
sius [^) that the rate of wall collisions per particle is 
Av/{'W), where A is here the area of the dividing plane, 
V is half of the containing volume and v is the speed of 
the particles, we can calculate the rate of passage between 
the two halves. From this we find 

tu, = —. (6) 

V 

Consider now a single UCN. Effectively, a coin is tossed 
once every to determine which side of the trap the 
neutron is in. Over a storage time t, this decision is 
therefore made N = t/ty^ times. The number of times n 
for which the UCN is on the side with the stronger field 
is binomially distributed with mean N/2 and variance 
7V/4. The additional J B ■ dt experienced by this UCN 
is ^ 2{n — N/2)ty;ABz (where the factor 2 accounts for 
the fact that when it is not in the stronger-field region 
it is in the weaker-field region). Multiplying this by the 
neutron gyromagnetic ratio 7 „ gives the extra precession 
angle Ot away from the mean. The polarization is the 
average projection upon the mean precession vector, and 
using 


+ (7) 

^2 

where T 2 is the transverse spin-relaxation time, together 
with COS0 ~ 1 — 0^/2 -I- ..., we find that 


Eq. in our calculations going forward, bearing in mind 
nonetheless that there is some uncertainty in the size of 
its contribution. The simulated results also appear in 
Fig. 2 of [T], where it is shown that the expected depen¬ 
dence upon {dBz/dz)^ holds true over a wide range of 
gradients. 



FIG. 2: (Color online) Intrinsic depolarization times T 2 for 
UCN in lOpT/cm magnetic-field gradients. The upper set 
are for vertical gradients dBzjdz with the results of simula¬ 
tions for completely diffuse (green circles) and 80% specular 
(green squares) reflections; the solid blue line is the analytical 
approximation of Eq. The lower set are for transverse gra¬ 
dients in the vertical field, of the form dB^/dx. Red upwards- 
pointing (downwards-pointing) triangles are for diffuse (80% 
specular) reflections, and the solid blue line represents the 
analytical approximation of Eq. |14[ T 2 scales as the inverse 
square of the applied gradient. 


After a measurement time t, the polarization a is re¬ 
duced by a factor 





( 10 ) 


t 




2,vgi 





( 8 ) 


where the subscript “vgi” stands for “vertical gradient, 
intrinsic”. This yields 


^lABlty, W^lidB^jdzr 

The upper solid blue line in Fig. shows the predic¬ 
tion of Eq. 1^ and despite the rather crude nature of its 
derivation it is seen to lie nicely between the results of 
our simulations for completely diffuse reflections (green 
circles) and for the case where the probability of specu¬ 
lar reflections is 80% (green squares). We therefore use 


which implies a parabolic profile to the dependence of a 
upon the vertical magnetic-field gradient. We see this 
in Fig. which shows (dotted black line) the spectrum- 
weighted average a as a function of the applied vertical 
gradient. 

We can now proceed to make a rough estimate of the 
extent to which this intrinsic depolarization may result 
in a shift in the measured neutron frequency. Intuitively, 
for example, it might seem that if low-energy UCN depo¬ 
larize more quickly than their high-energy counterparts, 
they would have less of a role to play in determining the 
frequency (since the uncertainty on the frequency mea¬ 
surement is inversely proportional to the polarization a). 
In consequence, the frequency measurement may appear 
to arise from a somewhat stiffer spectrum than is actually 
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shall see in Section VII| below) the Ramsey wrapping will 
in any case long since have taken hold. 


B. Horizontal gradients: dB^/dx 


We now carry out exactly the same calculations for 
the case of horizontal changes dB^/dx in the vertical 
magnetic field. The cylindrical trap is in this case to 
be bisected by a vertical rather than a horizontal plane, 
giving 


7rr 



(13) 


This then yields 


FIG. 3: (Color online) Polarization a, and fractional decrease 
in the shift 5R in frequency ratio, arising from intrinsic depo¬ 
larization only, after 180 s of storage. Black dotted (green dot- 
dashed): a from dBzjdz (dB^/dx). Blue solid (red dashed): 
relative decrease in SR, i.e. in the slope of R vs. dB^/dz, for 
vertical (horizontal) gradients in B^. Typically, dBzjdx may 
be 30-50 pT/cm. 


the case, thereby raising the effective center of mass of the 
neutrons and reducing the factor Ah in Eq. In order 
to consider this, let us for the time being imagine that we 
can make absolute frequency measurements without the 
complication of the modulo 2tt arising from the Ramsey 
measurement, thus ignoring the Ramsey wrapping that is 
characteristic of the behaviour of the lowest-energy UCN 

m- 

By definition, the height difference Ah between the 
centers of mass of mercury and UCN (Eq. is 


If zp{z,e) dz de 
JJ p{z, e) dz de 


( 11 ) 


where ZHg is the center of mass of the mercury atoms. 
We note in passing the standard result that, for e ^ H, 
Ah = i/^/(24e). When e < H, Ah = H/2 — QAe instead; 
in intermediate regimes. Ah may be derived from a more 
precise expression for the center of mass [Tl I14|. 

We now define an effective height difference 


. ff za{efp{z,e)dzde 


( 12 ) 


which takes into account the relative contribution of 
each energy bin to the frequency measurement. The 
intrinsic-depolarization induced fractional decrease in the 
frequency ratio, away from that anticipated by Eq. is 
then A{SR)/SR = {Ah — Ahes)/Ah. This function is 
shown (solid blue line) in Fig. We see that the effect 
stays at the 2% level or below until quite large vertical 
gradients, in excess of 500pT/cm, by which time (as we 


^2.hgi ~ 8r372(a^^/5a.)2’ (^ 4 ) 

where the subscript “hgi” indicates that this contribution 
arises from intrinsic depolarization due to the horizontal 
gradient. 

Once again, this result (lower solid blue line in Fig. 
provides a very resonable approximation to our 
simulations (red triangles - upwards-pointing, diffuse; 
downwards-pointing, 80% specular). The ^1.5 orders of 
magnitude difference in response between the vertical and 
horizontal gradients arises principally because the bottle 
is four times wider than it is tall, and the respective di¬ 
mensions enter to the third power. Fig. shows (green 
dot-dashed line) a as a function of the applied horizontal 
gradient. 

Within our apparatus, it is difficult to tune the hor¬ 
izontal gradient in Bq to better than about 30pT/cm, 
corresponding to 1 nT (one part per thousand) difference 
from one side of the bottle to the other. We note that 
50pT/cm would yield a T 2 of 700 s, perfectly consistent 
with that typically observed in the actual experiment and 
able to explain the reduction from an initial polarization 
of a = 0.86 when the trap is first filled to a = 0.67 at 
220 s storage time in the absence of a vertical gradient. 

We can also calculate, just as we did for the vertical 
gradient, the fractional decrease in the frequency ratio R 
that we might expect to see as a result of the spectral de¬ 
pendence of the intrinsic depolarization in this horizontal 
magnetic-field gradient. This is shown as a red dashed 
line in Fig. We see that at horizontal gradients of 
around 50pT/cm, the slope of R vs. dB^jdz (Eq. de¬ 
creases by about 10%. This factor would be a constant, 
independent of the applied vertical gradient - it would 
not impose any curvature upon the vertical-gradient de¬ 
pendence. However, these calculations are for illustration 
only: we remind the reader that the frequency averaging 
implicit here is invalid when using the Ramsey resonance 
technique. 
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C. Horizontal fields: dBt/dx 


We consider here additional weak fields Bt that are ev¬ 
erywhere parallel to the xy plane. If uniform, such fields 
simply act to produce a small tilt in the net direction of 
the main holding field Bq, and - since the perpendicu¬ 
lar components add quadratically - a tiny change in its 
magnitude. The resulting field would still be uniform, 
leaving both the depolarization rate and the frequency 
ratio R unaltered. 

Such horizontal fields may of course have gradients of 
their own, e.g. if they are quadrupole-type fields of the 
form Bx = qy, By = qx, as discussed in Section VI.c of 
[5] . The direction of the total field will alter slightly from 
one side of the bottle to the other, but the UCN spins 
follow these changes adiabatically during their trajectory. 

To understand the process in simple terms, let us first 
consider a neutron polarized with its spin along the z 
axis. If the cell has a difference ABt in a transverse 
(i.e. horizontal; x or y) field component from one side 
to the other, then on traversing the cell the UCN sees 
the B field tilt through an angle (j) = ABt/B^ in a time 
tc = 2rjvt, where Vt is the relevant transverse velocity. 
The angular frequency of this tilting motion is therefore 
Wtiit = {ABt/B^) ■ {vt/2r). To keep B steady, we go to a 
reference frame rotating at w/ = Wtut. To see the correct 
spin motion in this frame, we have to add the field 


_ 1 _ ABt Vt 

* B^ 27 „r' 


(15) 


A new B't must be used after each wall collision, since 
ojf changes abruptly at that point. The result is that 
the spin of any one UCN executes a random walk, trac¬ 
ing out cones of small opening angles 0 i, 02 ,^ 3 ---, where 
6 = B't/B^, in the vicinity of the z direction. Assuming 
N such wall collisions during a storage time t = Ntc, 
these small angles add vectorially to give a total angular 
displacement of magnitude 




2r B, 


(16) 


the collision frequency) we note that its dependence upon 
parameters is identical to that of the rather more sophis¬ 
ticated Eq. 66 in m- 

This particular depolarization mechanism is less effec¬ 
tive by a factor of two when acting upon a spin process¬ 
ing in the horizontal plane rather than aligned with z, 
for the simple reason that B^ components do not affect 
the X component of spin, and similarly for y components. 
This leads to the well known result for this case m 

T2,tfi = (18) 


Putting in realistic numbers for our apparatus (a few 
nT for ABt, and v ~ 2m/s), we find that Eq. 17 predicts 
Ti tfi (and therefore '72,tfi) values of order 10° seconds. 
We certainly cannot expect to be sensitive to this. In 
any case, we have observed Ti times in excess of 1000 
seconds, implying that T 2 ^tfi > 2000 seconds, so this is 
clearly not a dominant effect. 

In terms of frequency shifts, the mercury atoms will 
average out the horizontal components, whereas the neu¬ 
trons remain sensitive to the total field magnitude. This 
gives rise to a change in the frequency ratio R of 


SR q^r"^ 

"R ^ 


(19) 


where, as before, the radius of the trap is r. This will be 
a constant shift, independent of the applied vertical gra¬ 
dient. In the case of the EDM spectrometer, where hor¬ 
izontal field components are several hundred to a thou¬ 
sand times smaller than the vertical field, the resulting 
frequency shifts are of the order of a part per million or 
less. 


D. Wall collisions 

The cell walls may contain tiny magnetic impuri¬ 
ties. Collisions with these would disturb the spins on 
timescales much shorter than the Larmor precession pe¬ 
riod. Since such perturbations can affect any orientation 
of spin equally, one can anticipate that 


Following the same methodology as for Eq. and sub¬ 
stituting for B't from Eq. we arrive at 


Bt 

AS/ ’ 


(17) 


where we now refer to the longitudinal spin-relaxation 
time Ti rather than T 2 because we began with the spin 
aligned along z. The subscript “tfi” refers to “transverse- 
field, intrinsic”, and we have taken vj = u^/3. This 
derivation is of course extremely simplistic (for example, 
use of the mean free path A rather than the cell diame¬ 
ter would immediately reduce by a factor ~ 3 for 
our trap geometry). However, it gives interesting insight, 
and (given that we are in the “high-held” regime where 
the spin-precession frequency is substantially higher than 


T2,waU = Ti^wall- (20) 

As noted above, we have measured Ti to be in excess of 
1000 s, which therefore sets a lower limit of 1000 s on the 
contribution to T 2 arising from wall collisions. This is 
therefore unlikely to be a signihcant source of depolar¬ 
ization. 


V. POLARIZATION VS. APPLIED VERTICAL 
GRADIENT 

Having established the expected response to the intrin¬ 
sic depolarization, we now go on to look at the effects of 
the gravitationally enhanced depolarization, using calcu¬ 
lations as discussed in Section imi above. 
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We show in Fig.|^the residual polarization, after 180 s 
of storage, as a function of the applied vertical magnetic- 
field gradient. The data points (black triangles) represent 
measurements [5] made with the magnetic holding field 
Bq pointing downwards. These measurements were made 
in 2012, more than two years before the spin-echo mea¬ 
surements of the UCN spectrum, but we have no reason 
to suspect that the spectrum would have changed during 
the intervening period. The solid line in light magenta 
shows the approximate expected contribution of the in¬ 
trinsic depolarization, based on the formulae of Eq.|^and 
Eq.[^ its profile should be correct, although there is un¬ 
certainty over its scale because we do not know the extent 
to which reflection of UCN within the trap is specular. 
The blue solid (green dashed) line shows the contribu¬ 
tion of gravitationally enhanced depolarization, using the 
measured SE (MC) spectrum as input. The blue dotted 
(green dot-dashed) line shows the combined calculated 
effect. In each case the calculated profiles are normal¬ 
ized to reflect the peak measured value of 0.67, which is 
a result both of imperfect initial polarization and also of 
intrinsic depolarization e.g. from horizontal field gradi¬ 
ents of the form dBzjdx. Bq up data are omitted from 
this plot for clarity; they are similar in form to the data 
shown, but with a somewhat lower maximum value of 
about 55%. 



FIG. 4: (Color online) Depolarization as a function of applied 
vertical magnetic-field gradient. Light magenta line: calcu¬ 
lated intrinsic contribution. Blue solid and green dashed lines 
are the enhanced contributions applied to the SE and MC 
spectra, respectively. The blue dotted and green dot-dashed 
lines show the combined effects. Measured data points are 
represented as black triangles. 

There are two striking features about this plot. The 
first is the very distinctive peaked shape of the profile 
near the maximum. The intrinsic depolarization mecha¬ 
nism has a very soft peak, parabolic in nature. In con¬ 
trast, the gravitationally enhanced component is almost 


triangular in form, precisely mirroring the behaviour of 
the data. 

The second feature of interest is the close match in the 
polarization profile across a wide range of gradients. No 
parameters were optimized in the calculated curves be¬ 
yond the normalization of the peak value (equivalent to 
assuming dB^/dx = 50pT/cm). As noted above, the MC 
spectrum is expected to be a little too hard, and we see 
that the data lie below the corresponding (green) lines as 
one would anticipate. The measured SE spectrum, on the 
other hand, is known to be a little too soft, since it is rep¬ 
resentative of a 220 s storage time whereas the data points 
were measured at 180 s. One would therefore expect the 
combined effects of intrinsic and enhanced contributions 
(dotted blue) to lie a little below the data points, as in¬ 
deed they do. It would appear that, rather fortuitously, 
the offset from the use of a softened spectrum is here al¬ 
most exactly compensated by the additional contribution 
of the intrinsic depolarization, leaving the calculated en¬ 
hanced contribution more or less perfectly aligned with 
the data. 


VI. FREQUENCY SHIFTS AT LARGE 
VERTICAL FIELD GRADIENTS 

We now turn to measurements of the ratio of neutron 
to mercury precession frequencies under applied vertical 
magnetic-field gradients. 

Fig.|5] shows the measured data alongside the results 
from calculations of the effect of gravitationally enhanced 
depolarization, based upon the measured SE (solid blue 
line) and simulated MC (dashed green line) spectra dis¬ 
cussed in Section |IT] above. The adjacent blue dotted 
and green dot-dashed lines include the approximate re¬ 
spective contributions from intrinsic depolarization, with 
the sines and cosines of the contributing phases (see Eq. 

weighted as a^. As anticipated, the intrinsic depolar¬ 
ization has little additional effect. 

Both the SE and the MC spectra result in the right gen¬ 
eral trend, i.e. curvature of the appropriate form. How¬ 
ever, as one might expect, the (stiffer) simulated spec¬ 
trum results in curves that are less steep than the data, 
whereas the (softer) spin-echo spectrum, with its larger 
Ah, results in curves that are rather steeper. 

Also shown in this figure as a pair of red dashed lines 
is the expected response based on Eq. using the Ah 
from the SE spectrum (with no depolarization), as well 
as (black dot-dashed lines) a fit to the data points with 
vertical gradients of less than 60 pT/cm. Bearing in mind 
both that the SE spectrum is a little too soft and that 
no depolarization effects at all are included, one would 
anticipate that the former lines would be slightly steeper 
than the latter, as is indeed the case. 





FIG. 5: (Color online) Ratio of neutron to mercury fre¬ 
quencies as a function of the applied vertical magnetic-field 
gradient. Triangular upwards-pointing data points (running 
downwards diagonally) are for Bo up; triangular downwards- 
pointing data points (running upwards diagonally) are for Bo 
down. Red dashed lines show the expectation from Eq. 
using Ah from the SE spectrum. The black dot-dashed lines 
represent a fit to data in the central region, with gradients 
of less than 60 pT/cm. The blue solid (green dashed) line in¬ 
cludes gravitationally enhanced depolarization, based on the 
SE (MC) spectrum, with the adjacent blue dotted (green dot- 
dashed) lines including the effect of intrinsic depolarization. 



EIG. 6: (Color online) A closeup of the central region of 
the curves showing the dependence of the ratio of neutron to 
mercury frequencies upon the applied vertical magnetic-field 
gradient. Triangular upwards- and downwards-pointing data 
markers once again represent measurements taken with the 
holding field Bo aligned vertically upwards and downwards 
respectively. Solid blue and dashed green lines again arise 
from calculations based on the spin-echo and simulated spec¬ 
tra, respectively. The black dot-dashed lines in the corner 
regions represent the trend of the data at higher gradients. 


VII. FREQUENCY SHIFTS AT SMALL 
VERTICAL FIELD GRADIENTS 

The effect that we wish to discuss here is arguably 
more subtle. We focus upon the very central region of the 
frequency-ratio curves, as shown in Fig.|^ We again show 
(via the blue solid and green dashed lines, respectively) 
the results of calculations of gravitationally enhanced de¬ 
polarization based upon the SE and MC spectra. 

It is apparent that there is a change in the slope of 
the lines as the crossing point is approached. It is visible 
in the data, where the trend at higher gradients is high¬ 
lighted by the black dot-dashed lines: these represent a 
common fit of all of the data with gradients of more than 
lOOpT/cm to the function 


I i? — i?o I = m 


dB, 

dz 


+ c, 


( 21 ) 


where Rq = 3.8424574 is the crossing-point value of R 
[8]. If extrapolated to lower gradients, these lines would 
clearly result in a significant discontinuity. 

We ascribe this phenomenon to Ramsey wrapping, 
which makes the spreading low-energy tail of the array 
of integrated phases indistinguishable from the contribu¬ 
tions of high-energy UCN. This moderates the frequency 
shift as described in [T] and shown in Fig. 4 therein. Intu¬ 


itively, we would expect this to start happening when the 
difference in magnetic field between the bottom of the 
bottle (where the low-energy UCN preferentially spend 
their time) and the center of the bottle (which is the av¬ 
erage position for higher-energy neutrons) is sufficient to 
produce a phase shift of between tt and 27r over the 180 s 
storage time. This amounts to about 100-200 pT over the 
6 cm half-height, i.e. a gradient of about 15-30 pT/cm, 
which is precisely where we observe it happening. 

The particularly keen-eyed reader may be able to per¬ 
ceive a further slight curvature in the calculated curve at 
a gradient of about ±50 pT/cm. This is due to that same 
low-energy tail wrapping itself around for a second time. 
The data do not have adequate resolution to discern this 
effect. 

We note finally that both the curves and the real phys¬ 
ical behavior of the system are expected to be symmetric 
with respect to the crossing point (indeed, the calcula¬ 
tions represented here were generated in a single quad¬ 
rant only, and then reflected through dBzjdz = 0 and 
through R = Rq). Therefore, data taken at points using 
the same positive and negative gradients and with both 
Bq field directions would allow one to extract the correct 
crossing point without knowing the curvature. This is an 
important safeguard for future nEDM data taking. 








9 


VIII. IMPLICATIONS FOR NEDM 

The leading systematic error in the nEDM measure¬ 
ment described in arose from shifts in the Larmor pre¬ 
cession frequency brought about by the interplay between 
(a) small magnetic-field gradients within the apparatus 
and (b) the motional magnetic fields due to the particles 
(both UCN and, in particular, the mercury atoms used 
for magnetometry) moving through the applied electric 
field [3]. This effect was compensated by considering the 
behaviour of the observed apparent EDM signal over a 
range of magnetic-field gradients. There was no direct 
measurement of the field gradient: instead, it was pa¬ 
rameterized (Eq.[^ by the ratio R of neutron to mercury 
precession frequencies. It is now clear, however, that even 
at quite moderate gradients R is subject to the nonlin¬ 
earities discussed above. Furthermore, it is stated in |1] 
that the height difference A/i between the UCN and the 
mercury atoms was 2.8 mm, with a precision of 4%. This 
latter was based upon measurements of the frequency re¬ 
sponse within a variable-height trap; in fact those mea¬ 
surements also would have been affected by gravitation¬ 
ally enhanced depolarization. The UCN energy spectrum 
was undoubtedly softer than had been thought, implying 
that the calculated slope of the lines in Fig. 2 of |3] was 
steeper than it should have been. The relatively good 
match of the calculated to the actual slope of the fitted 
line is due in part to the nonlinear nature of the depen¬ 
dence of R upon dBz/dz, which would be similar in form 
to that shown in Fig. above and in Fig. 4 of [I], and 
which would result in a steeper-than-expected slope once 
beyond the linear region. 

Since the data were taken more or less symmetrically 
about the crossing point, this effect is unlikely to produce 
a very substantial change in the nEDM limit in this case. 
Nonetheless, a detailed reanalysis is now underway, and 
it is expected that a revised result will emerge shortly. 
Future measurements will doubtless enjoy the advantage 
of better diagnostics both of the magnetic field (with im¬ 
proved magnetometry) and of the energy spectrum (using 
the spin-echo technique). 

IX. CONCLUSIONS 

Measurements undertaken at the EDM spectrometer 
at PSI, showing the dependence upon applied vertical 


magnetic-field gradients of depolarization rates and of 
the neutron precession-frequency, have clearly demon¬ 
strated features that are characteristic of the anticipated 
behaviour resulting from gravitationally enhanced de¬ 
polarization and Ramsey wrapping: namely, a sharply- 
peaked rather than parabolic depolarization profile, and 
significant nonlinearities in the frequency-response curve. 
Using estimates of the spectrum of stored UCN, based 
upon measurements using the spin-echo technique and 
also upon detailed simulations, we have demonstrated 
excellent qualitative agreement between measurements 
and theoretical expectations. It also seems clear that 
intrinsic depolarization processes have only a marginal 
effect upon frequency shifts in the presence of magnetic- 
field gradients, and that such shifts are dominated by the 
gravitationally enhanced component. 

There are obvious implications for nEDM measure¬ 
ments, including for the analysis that led to the current 
world limit |1] , since the frequency-response curve is used 
to correct for systematic effects. 
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